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We study coherent backscattering phenomena from single and multiple stacking faults (SFs) in 
3C- and 4H-SiC within density functional theory quantum transport calculations. We show that 
SFs give rise to highly dispersive bands within both the valance and conduction bands that can 
be distinguished for their enhanced density of states at particular wave number subspaces. The 
consequent localized perturbation potential significantly scatters the propagating electron waves 
and strongly increases the resistance for n-doped systems. We argue that resonant scattering from 
SFs should be one of the principal degrading mechanisms for device operation in silicon carbide. 
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Silicon carbide is a wide band-gap semiconductor with 
a formidable number of stable polytypes (> 200) that are 
distinguished by a unique stacking sequence of SiC bi- 
layers. The origin of this vast polymorphism stems from 
the small differences in the formation energies between 
different stacking sequences that outscore all other semi- 
conductors by at least one order of magnitude^! Likely, 
the same reason lies underneath the frequent presence of 
stacking faults (SFs), i.e. irregularities in the stacking 
sequence, found in this material. Even if the structural 
lattice deformation induced by SFs is marginal, their 
electronic response is that of localized defects; they in- 
troduce particular quantum-well-like states in the eigen- 
spectrum, with wave functions that are confined within 
~10 A around the SpEHH. Since the early days of SiC 
device fabrication, there was a clear evidence that such 
defects were related to a degradation of the device op- 
eration, usually giving leakage currents during electrical 
conduction^. However, such vision overlooks the possibil- 
ity of SFs acting as scattering centers during the trans- 
port process^. From a fundamental point of view, the 
presence of a localized potential barrier in a material (e.g. 
from a defect or impurity) can either partially or totally 
reflect the incident electronic waves, on the basis of its 
height and energy resonance^. Stacking faults can be 
no exception to this. Rather, their plausible microme- 
ter dimensions can make this aspect crucial. In this case 
the perturbation induced on the transmitting waves is 
not localized at an atomic position (like in the case of 
vacancies, interstitials, etc.), but it is spatially extended 
throughout the entire width of the dislocation. 

In this article we study the coherent backscatter- 
ing mechanism induced by stacking faults in cubic and 
hexagonal SiC within a computational approach that 
treats the electronic structure from first principles (us- 
ing the density functional theory) and the quantum 
transport within the nonequilibrium Green's function 
formalism^]. We show that the effect of this scat- 
tering mechanism can be extremely important for the 
conduction-band transport characteristics, giving rise to 
total electronic reflections for particular resonances and 
types of defects. 
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FIG. 1. Crystal structure of 3C-SiC in the presence of SFs: 
(a) ideal, (b) single SF, (c) double SF and (d) triple SF (micro- 
twin) . Transverse (e) and parallel (f ) views with respect to the 
transport direction, showing the lattice vectors of the hexag- 
onal unit cell. 



We commence by simulating three-dimensional 3C-SiC 
systems with single, double and triple (micro-twin) SFs 
perpendicular to the [111] crystallog raph ic direction [Fig. 
[IJa-d)] using the supercell methooP^. We construct 
the supercells using a hexagonal symmetry for the (111) 
plane while the c-axis is parallel to the [111] direction 
[Fig. |lje-f )] . The c-axis also defines the direction of 
transport, making the incident electron angle <fi with re- 
spect to the SF equal to <f> = 90° . We consider a single SF 
within the supercell, ensuring that the cell length guar- 
antees the absence of SF interactions between successive 
periodic images. We calculate the electronic structure us- 
ing the SIESTA code^ along with the Local Density Ap- 
proximation (LDA]f31 for the electronic correlations, in 
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order to have a clear comparative framework with Ref. 
I -5 1 G 1 1 12 L The wave functions are written on a localized 
basis set of double-^ valence plus polarization atomic or- 
bitals for both Si and C, while the electronic contribu- 
tion of the ionic cores is described with norm-conserving 
Troulicr-Martins pseudopotentials^. The minimization 
of the electron density is achieved by sampling the Bril- 
louin zone with an 8 x 8 x 1 Monkhorst-Pack grid, whereas 
a mesh cutoff energy of 450 Ryd is used for real-space in- 
tegrals. Structural relaxation is obtained with a force 
criterion of 0.04 eV/A. Similar calculations are also con- 
ducted for silicon. 

Upon convergence, the Hamiltonian (H) and Overlap 
(S) matrices are extrapolated for 140 nonequivalent k- 
points along the T —> M —> K — > T Brillouin-zone path, 
which is parallel to the SF junction (fcy from now on). 
Quantum transport is calculated for the ballistic regime, 
starting from the single-particle retarded Green's func- 
tion matrix 

G r h (e) = [eS k[l -H h - E £>fc „ - S^J" 1 , (1) 

where e is the energy and are self-energies that 

account for the effect of ideal semi-infinite contacts (ob- 
tained from surface Green function techniques^-). The 
transmission probability of an incident Bloch state with 
energy e and momentum fen can be thereon computed as: 

T fe|| (e)=Tr-{r L , fe|| ^ || r fl , fe|| K || ] t }, (2) 

where T L ( R ) ik]] = i{£L(R),k {l ~ p£(fl),fe|| ] f } are the 
spectral functions of the contacts. Consequently, the re- 
flection coefficient can be defined as R = 1—T. According 
to the Landauer-Buttiker theory^, conductance can be 
calculated as G = GqT, where Go = 2e 2 /h « 77.5/iS is 
the conductance quantum. Similarly, from the spectral 
function of the scattering region Afc,, = i\Q\. — ]t} 
the k-resolved density of states (DOS) can be computed 
as: 



DOS h (e) = Tr[A h S h /(2*)] (3) 

It has to be noted here that the total transmission and 
DOS in real space can be estimated for each energy by 
integration over the reciprocal space. In this case, the 
number of the sampled fc-points for each Brillouin zone 
subpath becomes proportional to the dimension of the 
corresponding real-space direction, being a measure of 
confinement for the crystal. 

We start by considering a momentum-space projection 
of the main transport attributes of 3C-SiC. The elon- 
gated structure of the supercell along the [111] direction 
gives rise to a quasi-two-dimensional Brillouin zone with 
the r point lying at its center, whereas the M and K 
points are positioned at the center and the angle of the 




FIG. 2. Momentum-space projection of the conductance as a 
function of energy (in units of Go = 2e 2 /h) for (a) ideal 3C- 
SiC, (b) 3C-SiC with a single SF, (c) 3C-SiC with a double 
SF, (d) 3C-SiC with a micro-twin SF, (e) ideal Si and (f) Si 
with a singe SF. 



hexagon's side respectively^. The fc-resolved ideal con- 
ductance of pure 3C-SiC [Fig. |2[a)] has some partic- 
ular characteristics with respect to these high-symmetry 
points: (i) the valence band maximum is mapped at the T 
point whereas the conduction band minimum is mapped 
at the M point. Therefore these are the points of interest 
for low-bias p- and n-type conduction respectively, (ii) 
There exists only a single conduction channel per fc-point 
for the conduction band up to ^2 eV from its minimum, 
as can be deduced by the = Go value of the conduc- 
tance for this energy range. Contrary, the conductance 
value at the maximum of the valence band is 3Go- This 
feature makes cubic SiC particular within the SiC poly- 
types, since the breaking of the cubic symmetry results in 
higher values of the conductance at the minimum of the 
conduction band (e.g. in the case of 4H-SiC G/.,, = 2G 
near the conduction band minimum [see Fig. |4^a)]). The 
presence of a single SF [Fig. [2](b)] drastically changes this 
picture. Band structure calculations already exclude the 
possibility of intra-gap states from SFs in 3C-SiC, con- 
trary to other hexagonal polytypeiP. SF-induced states 
in this case are introduced inside the valence and con- 
duction bands. Their influence on the system is pertur- 
bative and as a result, the transport channels get par- 
tially blocked. The local effective potential induced by 
the single-SF scatters the electron waves mainly around 
the M point at the conduction band and around the K 
point in the valence band. Hence, scattering by the SF 
has an important effect on the low-bias conductance of 
n-type systems and a marginal influence on the conduc- 
tance of p-type systems. The percentage of reflection is 
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86% near the conduction band minimum, meaning that 
the single-SF barrier strongly increases the resistance, 
albeit not completely suppressing it. This picture quan- 
titatively changes in the case of the double and the triple 
SFs [Fig. |2jc-d)]. Here the increase of the barrier height 
due to the spatial extension of the defect in real space 
gives rise to an almost complete reflection of the electron 
waves at the bottom of the conduction band [see table 
[i], and transport effectively takes place only for higher 
energies. In this case, this extreme enhancement of the 
resistance can be directly related to a respective increase 
of the measured Schottky-barrier height for n-type 3C- 
SiC systems; indeed, even if the difference between the 
electronic states at the conduction band minimum and 
the Fermi level is not altered by the SF, electrical mea- 
surements are expected to measure the effective height 
between the Fermi level and the first states within the 
conduction band that do not totally reflect the electrons. 
Our calculations indicate that the theoretical limit for 
this increase could be as high as ~ 0.8eV with respect to 
the ideal case for charges approaching the SFs at right 
angles; however, considering the presence of electron de- 
phasing in real systems, the measured value is expected 
to be lower. For all the previous cases, transport at the 
top of the valence band practically remains unaffected, 
implying that resonant backscattering from SFs should 
not significantly influence conduction for p-type 3C-SiC 
systems. 

The origin of this strong scattering mechanism in 3C- 
SiC is inherently related to the geometrical symmetry of 
the SFs (being perpendicular to the [111] direction) with 
respect to the directionality of the atomic orbitals that 
form the wave functions at the minimum of the conduc- 
tion band (considering a linear wave function expansion 
on an atomic orbital basis). A direct comparison with 
the case of Si (ideal and with a single SF [Fig. [2^e-f )] ) 
is illuminating: a first reading shows that the main dif- 
ference between the two systems lies in their conduction 
bands, which are clearly inequivalent. If we now exam- 
ine the single-orbital contributions in the wave functions 
for each system^, we find that for Si, the conduction 
band minimum has principal contributions from atomic 
orbitals that are not polarized towards the [111] direc- 
tion. Hence, the SF does not break the symmetry of the 
respective wave function and as a result reflection is al- 
most absent. Contrary, in the case of 3C-SiC, orbitals 
with a directionality parallel to the c-axis do contribute 
in the formation of the wave functions at the conduction 
band minimum. In this case the geometrical symmetry 
breaking is reflected into an electronic symmetry break- 
ing, which forms the origin of the conductance pertur- 
bation. Moreover, for higher energies, silicon has numer- 
ous bands throughout the entire Brillouin zone, which all 
contribute in the total conductance of the system. Based 
on the polarization criterion, the states introduced by 
the SF affect only some of them, whereas others remain 
unaffected. As a result reflection is overall moderate. 

A key point that emerges with respect to the trans- 



System 


Reflection coeff. (%) 


E M (eV) 


3C-SiC single SF 


86 


0.1/1.9 


3C-SiC double SF 


98 


0.1/1.2 


3C-S1C triple SF 


100 


0.2/0.8/2.2 



TABLE I. (Central column) Reflection coefficient at the bot- 
tom of the conduction band. (Right column) SF energy res- 
onances from the bottom of the conduction band at the M 
Brillouin- zone point. 
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FIG. 3. Momentum-space projection of the density of states 
as a function of energy (in arbitrary units) for (a) 3C-SiC with 
a double SF and (d) 3C-SiC with a micro-twin SF. Arrows 
indicate the bands that originate from the SFs. 



port features presented in the previous paragraphs is the 
nature of the SF states in 3C-SiC. Previous studies^' 
have indicated that these states should have resonances 
at the bottom of the conduction band and excluded any 
intra-gap contribution. However their degree of local- 
ization and mixing with the non-localized bands of the 
crystal make their exact identification unclear for band 
structure calculations, and similarly difficult for exper- 
imental measurements. In our approach, based on an 
exact estimation of the density of states in the recip- 
rocal space (see eq. [3]), we show that SFs give rise to 
highly dispersive bands within both the conduction and 
the valence zones. These bands are characterized by an 
enhanced density of states within particular wave num- 
ber subspaces, which correspond to the Brillouin-zone re- 
gions of enhanced backscattering during conduction [Fig. 
|3^a-b)]. Table [I] shows the resonant energies of the var- 
ious SFs studied in the previous paragraph for 3C-SiC 
at the M Brillouin-zone point. A common aspect for all 
SFs is that they generate a state that is very close to the 
conduction band border with a "faint" DOS signature. 
Moving onwards to higher energies, SF bands are char- 
acteristic of the SF type and have a higher DOS with 
respect to the non-defected crystal. 

Results regarding cubic SiC show that the physical ex- 
tend of the stacking fault can be extremely important 
for its scattering properties, giving rise to electronic re- 
flections that range from partial to total. Qualitatively, 
the same mechanism lies underneath the coherent scat- 
tering process in hexagonal SiC polytypes. As an exam- 
ple, here we simulate the extended (4,4) Shockley-type 
defect!^ [Fig. Qc)] that appears to be one of the most 
common SFs in 4H-SiC crystals^. The methodology 
remains identical as in the case of 3C-SiC, considering 
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FIG. 4. Momentum-space projection of the conductance as a 
function of energy (in units of Go = 2e 2 /h) for (a) ideal 4H- 
SiC, (b) 4H-SiC with a (4,4) SF. (c) Schematic representation 
of the 4H-SiC supercell, highlighting the (4,4) SF. 



now that the SF is perpendicular to the [0001] direction. 
The computational outcome [Fig. |2]jb)] indicates a se- 
vere transport damage for the conduction band and a 
respective increase of the material's resistivity along the 
c-axis for n-doped systems. S uch fe ature is in excellent 
agreement with the experiment 19 l 2 ° l, indicating that res- 
onant scattering effects are indeed at the origin of the 
transport anisotropy in SiC crystals. Similar phenomena 
have been theoretically predicted for III-V nitride semi- 
conductors within a semiclassical picture^. 

A final issue that needs discussion is the weight of co- 
herent backscattering with respect to phase incoherency 
and other diffusive scattering mechanisms that exist in 
the electrical conduction of realistic 3C-SiC crystals. We 
can expect that dephasing of the electronic waves should 
flatten the conductance distribution and smear its max- 
ima and minima. Moreover, dissipative events (e.g. 
electron-phonon interactions) are expected to further in- 
crease the resistance with respect to the ideal case, tuning 
the electrical conduction according to Ohm's law. How- 
ever, the impact of the coherent backscattering mecha- 



nism for energies around the minimum of the conduction 
band is expected to be prevailing for extended SFs. The 
angle of the incident electron waves with respect to the 
SF is again an important factor for this phenomenon. 
The results shown here (cf> — 90°) represent a condition 
of maximum backscattering. For 90° < 4> < 0° gradual 
leakage currents are expected to pass through the SF bar- 
riers and increase the conductance. In the case of <t> = 0°, 
.i.e. when SFs are parallel to the transport direction, we 
can expect an absence of any resonant scattering effect. 
In this case the increased DOS of SFs should give rise to 
a higher localized current density with respect to that of 
the rest of the crystal (in the case when the SF is con- 
tacted by a metallic electrode) . Such characteristic could 
explain recent current maps obtained by n-type 3C-SiC 
structures using conductive atomic force microscopy^]. 

To conclude, in this letter we have investigated the 
electronic conduction properties of SFs in 3C- and 4H- 
SiC within first-principles quantum transport calcula- 
tions. We have found that SFs strongly scatter the elec- 
tron waves at particular energies and Brillouin zone sub- 
spaces for both the conduction and the valence band. 
However, a significant increase in the material's resistiv- 
ity is only expected for n-type systems and incident elec- 
tron angles other than (j> — 0° . In addition, the presence 
of extended SFs should totally reflect the electron waves 
at the bottom of the conduction band, effectively increas- 
ing the height of the measured Schottky-barrier. Based 
on these results we argue that, contrary to silicon, elec- 
tron resonant scattering from SFs should be one of the 
principal sources of instability and device degradation in 
silicon carbide. 
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